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Establishing the neutrino mass hierarchy is one of the fundamental questions that will have to be 
addressed in the next future. Its determination could be obtained with long-baseline experiments 
but typically suffers from degeneracies with other neutrino parameters. We consider here the NO^A 
experiment configuration and propose to place a second off-axis detector, with a shorter baseline, 
such that, by exploiting matter effects, the type of neutrino mass hierarchy could be determined with 
only the neutrino run. We show that the determination of this parameter is free of degeneracies, 
provided the ratio L/E, where L the baseline and E is the neutrino energy, is the same for both 
detectors. 

PACS numbers: 14.60.Pq 



I. INTRODUCTION 

During the last several years the progress in the studies of neutrino oscillations has been remarkable. The experi- 
ments with solar pj 0, y> il- M- Q > atmospheric Q , reactor || and recently also long-baseline accelerator [!j neutrinos 
have provided compelling evidence for the existence of neutrino oscillations driven by non-zero neutrino masses and 
neutrino mixing. We know that there are two large (#12 and #23) and one small (#13) angles, and at least two mass 
square differences , Am^ = m 2 — mf , with m^i the neutrino masses, one associated to atmospheric neutrino oscilla- 
tions (Amjj) and one to solar ones (Am§i). The angles 6*12 and 023 represent the neutrino mixing angles responsible 
for the solar and the dominant atmospheric neutrino oscillations, while $13 is the angle limited by the data from the 
CHOOZ and Palo Verde experiments |l3l Il4f . 

Stronger evidences of neutrino oscillations are provided by the new Super-Kamiokande data on the L/E dependence 
of multi-GeV /i-like atmospheric neutrino events |15|. L being the distance traveled by neutrinos and E the neutrino 
energy, and by the new more precise spectrum data of the KamLAND ^(| and K2K experiments 9] . For the first 
time these data exhibit directly, not only a deficit with respect to the expected si gnal , but also the effects of the 
oscillatory dependence on L/E of the probabilities of neutrino oscillations in vacuum We begin to actually "see" 
the oscillatory behavior of neutrino propagation. 

The Super-Kamiokande and K2K data are best described in terms of dominant — > v T [p^ — > v T ) vacuum 
oscillations. The best-fit values explaining the Super-Kamiokande data Q are lAm^l = 2.1 x 1CP 3 eV 2 , sin 2 26*23 = 
1.0 , whereas those for the K2K data @ are |Am| x | = 2.8 x 10" 3 eV 2 , sin 2 26* 23 = 1.0 . The 90% C.L. allowed ranges 
of the atmospheric neutrino oscillation parameters obtained by the Super-Kamiokande experiment read : 

lAm^l = (1.5-3.4) x 10~ 3 eV 2 , sin 2 26> 23 > 0.92. (1.1) 

The sign of Am| x and of cos 2023, when sin 2 26*23 7^ 1-0, cannot be determined with the existing data. For the 
mass square difference, the two possibilities, Ato 2 i > or Am^j < 0, correspond to two different types of neutrino 
mass ordering: normal hierarchy (NH), mi < ui2 < TO3 (Amf^ > 0), and inverted hierarchy (IH), 7713 < mi < TO2 
(Am| x < 0). The fact that the sign of 008 26*23 is not determined when sin 2 26*23 7^ 1-0 implies that the octant where 
6*23 lies is not known. 

In addition, the combined 2-ncutrino oscillation analysis of the solar neutrino data, including the results from the 
complete salt phase of the Sudbury Neutrino Observatory (SNO) experiment [||, and the recent KamLAND 766.3 



1 We restrict ourselves to a three-family neutrino scenario analysis. The unconfirmed LSND signal 111! cannot be explained in terms 
of neutrino oscillations within this scenario, but mi ght require additional light sterile neutrinos or more exotic explan atio ns (see e.g. 
Ref. Illl ). The ongoing MiniBooNE experiment 11 J is expected to explore all of the LSND oscillation parameter space I 111 . 
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ton-year spectrum data shows that the solar neutrino oscillation parameters lie in the low-LMA (Large Mixing 
Angle) region, with the best fit value at @ 

Amli = 8.0 x 1CT 5 eV 2 , sin 2 12 = 0.31. (1.2) 

In a 3-neutrino oscillation framework, a combined analysis of the solar, atmospheric, reactor and long-baseline neutrino 
data gives (see also Ref. 0]): 

sin 2 6»i3 < 0.041, 3cr C.L. (1.3) 

As we know that neutrinos do oscillate, there are very important questions that will have to be addressed in future 
experiments. Besides the more accurate determination of the leading neutrino oscillation parameters that will be 
achieved by the MINOS [20|]. OPERA [2lJ and ICARUS [2^ experiments and future atmospheric and solar neutrino 
detectors 2j|, one of the most important tasks in the next future will be the determination of the (1,3) sector 
of the lepton mixing matrix, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) neutrino matrix 24] . A complete 
determination of this sector entails the measurement of a non-zero #13, which will open the door to the experimental 
measurement of the CP- (or T-) violating phase, S, and to establishing the type of neutrino mass spectrum. This 
mixing angle controls the — > v e and — > T> e conversions in long-baseline experiments and in the widely discussed 
very long-baseline neutrino oscillation experiments at neutrino factories |2^,|2^|. More recently, /3-beams experiments, 
exploiting neutrinos from boosted-ion decays |27l l28l l29j | , with an improved experimental setup have been shown 
to achieve sensitivities to leptonic CP-violation and to the sign of the atmospheric mass difference competitive 
with those at neutrino factories |3(j. The mixing angle #13 also controls the Earth matter effects in multi-GeV 
atmospheric [ail El HI E3 

and in supernova 36] neutrino oscillations. Finally, the magnitude of the T-violating 
and CP-violating terms in neutrino oscillation probabilities is directly proportional to sin #13 j3?|]. Therefore the 
determination of the magnitude of #13 is crucial for the future searches for matter effects and CP-violation in the 
lepton sector at neutrino oscillation experiments. 

The measurement of, or a stron ger limit on , th e mixing angle #13 in the near future is going to be achieved by 
reactor [3^ and long-baseline pol l2ll I22I I3H. |40L l4l| neutrino experiments. Neutrino reactor experiments, being 
disappearance experiments, are only sensitive to the value of #13. Long-baseline neutrino experiments, in addition to 
having a better sensitivity to #13, are the only way (in the near future) to search for CP violation, while being able 
to determine the type of hierarchy at the same time 2 . However, long-baseline neutrino oscillation experiments suffer 
from degeneracies in the neutrino parameter space |50l IFll I52I I53I l54| . In general, the proposed experiments have a 
single detector with the beam running in two different modes, neutrinos and antineutrinos. With only one neutrino 
and one antineutrino run, the degeneracies can lead to different CP-violating and CP-conserving sets of parameters 
explaining the data at the same confidence level. In Ref. j2|| it was pointed out that some of the degeneracies could 
be eliminated with sufficient energy or baseline spectral information. In practice, however, the spectral information 
has been shown to be not strong enough to resolve degeneracies with a single detector, once that statistical errors and 
realistic efficiencies and backgrounds are taken into account. In order to resolve the par ameter deg eneracy, another 
detector [m1 |5^, |5(| 1^3, or the combination with another experiment [5^, IgJl EOl 

(H ilil El would, thus, be 

necessary. Recently, new approaches for determining the type of hierarchy have been proposed [66]; they typically 
exploit other neutrino oscillations channels, such as muon neutrino disappearance, and require very precise neutrino 
oscillation measurements. 

Contrary to the naive expectation, it has been shown numerically in Ref. |60| and analytically in Ref. |6l| that the 
use of only a neutrino beam could help in resolving the type of hierarchy when two different long-baseline experiments 
are combined under certain conditions. 

Differently from this approach, we present here a scenario with only one experiment, which runs in the neutrino 
mode and uses two detectors at different distances and different off-axis angles. It is well known that off-axis neutrino 
beams have a very narrow neutrino spectra, and that the peak energy can be tuned by just moving the detector with 
respect to the main beam axis. We notice that an off-axis beam can be obtained by either displacing the detector a 
few km away from the location of an on-axis surface detector, or by placing it on the vertical of the beam-line but at 
a much shorter distance. In such a way, a single beam could do the job of two beams with different energies. 



The same capabilities with better sensitivity could be achieved by neutrino factorie s |2fl |2fil . Information on the type of neutrino 
mass hierarchy might be obtained in future atmospheric neutrino experiments l32l l33t l33l3aT lf neutrin o ar e Majorana par ticles, next 
gen eration neutrinoless double j3— decay experiments could establish the type of neutrino mass spectrum |42| (see also Refs. |4 1 l44l l"tH 
|4(I ) an d, p ossibly, might p rovide some in formation on the presence of CP— violation in the lepton sector due to Majorana CP-violating 
phases |t3 (see also Refs. |3 EE El E3 ) • 
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We study the use of two off-axis detectors in combination with the NuMI beam so that the type of mass spectrum 
could be determined free of other degeneracies, if #13 is not very small. We will consider, for one of them, the location 
which is the most likely for the NOz^A configuration 3 (L — 810 km and E = 2.3 GeV) and show different possibilities 
for the baseline of the other detector in order to make the measurement of sig^Am^) feasible with only the neutrino 
beam. We name this improved experimental setup Super-NO^A [(57| • Following the line of thought of Ref. [§l|, we 
will show that a configuration with the same vacuum oscillation phase, i. e. same L/E for both detectors, is specially 
sensitive to matter effects. For such an experimental setup the sensitivity to sig^Am^) would be enhanced as the 
difference in baseline lengths grows. This configuration also has the advantage of requiring only one experiment and of 
reducing the error due to systematic uncertainties from the beam. In addition, we will show that such a measurement 
is free of degeneracies, which provides the full power of this method. We start by presenting the general formalism 
in Sec. [H] In Sec. II I II we describe the experimental setup. We explain in Sec. IIVI the method to extract the type 
of neutrino mass spectrum free of degeneracies by using this special configuration and we show how the sensitivity 
changes for different values of IAtttJ-jJ. Finally, in Sec.0 we make our final remarks. In Appendix 1X1 we present the 
computed charged-current (CC) neutrino event rates for a particular choice of parameters. 



II. FORMALISM 

We consider the probability of'j/'^ — > S? e oscillation, P fas' ^ -^V^L^ = P{L) 7 in the context of three-neutrino 

mixing. For neutrino energies E > 1 GeV, 6*13 within the present bounds [TsUli^. and baselines L < 1000 km |5^| 4 , 

the oscillation probability P (L) can be safely approximated by expanding in the sm all p arameters #13, A12/A13, 
A 12 /A and A 12 L , where A 12 = Am| 1 /(2£) and A i3 = Am^ 1 /(2E) [2fj (see also Ref. [68]): 

P\ L ) * sin 2 23 sin 2 20 13 (^) ' sin 2 (Mi^L ) 

+ cos 13 sin 20 13 sin 2023 sin 20i2 sin (4^) sin ( (i ^ 13)L ) cos (A^lL t S ) ( 2A ) 

+ cos 2 # 23 sin 2 26 12 ( Af ) 2 sin 2 (4p) , 

where L is the baseline. We use the constant density approximation for the index of refraction in matter A, defined 
as A = \/2GFn e {L), with n e (L) the average electron number density, defined by n e {L) = l/L J Q L n e (L')dL' . Here 
n e {L) is the electron number density along the baseline. 

As is well known [H^ . the CP trajectory in bi-probability plots of neutrino and antineutrino conversion at the same 
baseline, is elliptic under the assumption of mass hierarchy and of adiabaticity. The ellipses obtained for each of the 
two hierarchies, and for different values of 8, $13 and of 023, intersect in points that correspond to 2-, 4-, and 8-fold 
degeneracies |5fl l5ll l52l I53L l54| . It follows that even a precise determination of a point in the P-P plane can result 
in different sets of CP-conserving and CP-violating parameters (8, 0i3, 023, si(7n(Arn 2 1 )), all of which reproduce 
the observations. The allowed regions in the P-P plane obtained by varying the values of 0i3 and 8 within their 
ranges describe wide "pencils" . The "pencils" for the cases of normal and inverted hierarchy have different slopes 
and overlap for a large fraction. This indicates that, generically, a measurement of the probability of conversion for 
neutrinos and antineutrinos cannot uniquely determine the type of hierarchy in a single experiment. In order to resolve 
such degeneracy, various strategies have been proposed: combined analysis of the data from different super-beam 
experiments, or of the data from super-beam facilities and neutrino factories (or /3-beams) |30l l59l l6fl l6lt l62l l63| . 
the use of additional information from atmospheric neutrino data |64|. and experimental setups with clusters of 

detectors [mm mm. 

It has been pointed out that considering the probabilities of neutrino oscillations only, at two different baselines and 
energies, can resolve the type of hierarchy |6Ctl6l|. In the case of bi-probability plots of neutrino-neutrino conversions 
at different baselines, the CP-trajectory is elliptic too. Again, the allowed regions in these bi-probability plots form 
two "pencils" , which grow in width away from the origin, each of them associated with one type of mass spectrum. 
The overlap of the two "pencils" , which signals the presence of a degeneracy of the type of hierarchy with other 
parameters, is controlled by the slope and the width of the "pencils" . From Eq. (|2.1|l one can see that the ratio of 
the slopes of the central axes of these two "pencils" in the Pp-P^ plane, where Pp (Pn) is the neutrino conversion 



3 This values correspond to the already old NOi>A proposal |4(I . For the recent revised proposal see Ref. 1-411 . 

4 For E 0.6 GeV we have checked that the analytical expansion is accurate for L < 500 km within the present bounds of 6*13 and 
Amli/Amfi. 
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probability at the far (near) detector, is given by [6l| 

/ Aim \ 2 . 2 / (A n -A 13 ,n)-Ln N \ / _Ai^p \ 2 .9 ( (A f +Ai 3 , f )£f 
\A N -A 13 , N ) Sm ^ 2 ) \A F + A 13 , F ) Sm ^ 2 



a+ 



a (^|^) 2 S in 2 ( (^-^ ) (^t^) 2 sin 2 



(2.2) 



where a+ and a_ are the slopes for normal and inverted hierarchy, respectively; A 13 f(n): ^f(n) and £f(n) are the 
values of A13, A and L for the far (near) detector. Note that although we are using the constant density approximation, 
^4f and A-^ are different because the average density depends on the baseline. Using the fact that matter effects are 
small (A <§C A13), we can perform a perturbative expansion, which up to first order gives this ratio of slopes as 

a_ " \(Ai3,nLn/2) tan(Ai 3 ,N Ln/2)) 2AfLf ((Ai 3 ,fLf/2) tan(Ai3, F L F /2) 

For the case of L/E constant, noting that l/x — l/tanx is a monotonically increasing function, we conclude that 
the smaller the chosen energy, the larger the ratio of slopes. This ratio increases also for certain configurations with 
different L/E ,61). Another very important feature is the width, which is very small for equal L/E, but grows 
rapidly when this is not the case |6l| . Hence, even when the separation between the central axes of the two regions 
is substantial, unless the ratio L/E is kept close to constant, the ellipses overlap, making the discrimination of 
sig^Am^) challenging. As was shown in Ref. away from the L/i5-constant case, the choice L-p/E-p > L^/E^ 
is preferred. Otherwise, no matter how accurate the experiment is, the discrimination of the type of neutrino mass 
hierarchy free of degeneracies will not be possible. As a matter of fact, this is precisely what will happen when 
combining NOi^A and T2K experiments, for which -Lno^a/ -E^no^a = 352km/GeV < 421km/GeV =£t2k/£t2K 5 - In 
this case, a joint analysis of these two experiments does not present any interesting synergy effects and just accounts 
for adding in statistics 60]. 

Thus, we will consider the case of L/E constant and show how, by adding another detector to the already proposed 
NOi^A experimental setup, the measurement of the sign^Am^) is possible free of degeneracies from other parameters. 



III. EXPERIMENTAL SETUP 



As was pointed out in the previous section, we consider here only one experiment by using the same beam but 
having two detectors. In order to maximize the sensitivity to the type of hierarchy, we propose a configuration in 
which the neutrino conversion takes place predominantly at the same L/E at the two locations. In such a way, matter 
effects can be factored out and the type of neutrino mass hierarchy can be determined (if $13 is large enough), free of 
degeneracies from other parameters of the neutrino mixing matrix. 

In order to have the same L/E for both detectors we would need a very well peaked spectrum at both sites. This 
can be achieved by placing the detectors off the central axis of the beam. As is well known, most of the neutrinos in a 
conventional neutrino beam are produced in two-body decays tt^ — » /i* + < ^ > M . The energy and flux of these neutrinos 
is determined by the decay angle 9 |69|: 



where 7 is the Lorentz factor of the pion, 8 is the angle between the pion and the neutrino directions, and L is the 
distance between the decay point and the detector. A neutrino beam with narrow energy spectrum can be produced 
by placing the detector off-axis, i. e., at some angle with respect to the forward direction flbeam = 0. By using off-axis 
beams, one manages a kinematic suppression of the high energy neutrino components, whereas the low energy flux 
is kept approximately the same as that of the on-axis beams. The neutrino spectrum is very narrow in energy and 
peaked at lower energies with respect to the on-axis one. The suppression of the high-energy tail of the spectrum 



5 This corresponds to the configuration with an off-axis angle of 2° (OA2°) for T2K 39] . 
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greatly reduces the backgrounds due to neutral-current interactions and r production. Since the neutrino flux is 
nearly monochromatic, the off-axis technique allows a discrimination between the peaked v e oscillation signal and the 
intrinsic v & background which has a broad energy spectrum. An efficient reduction of the intrinsic background can 
therefore be achieved |69). 

The off-axis angle and baseline for each detector must be chosen in such a way that we have the same L/E at both 
sites. From Eqs. (|3 . 1|> and l|3.2|l we see that the flux scales as $ ~ (E/L) 2 ; so for this special configuration the flux 
at the near detector is of the same order as that at the far detector (see Appendix [SJ. For the far detector we will 
use the configuration proposed for the NOz/A experiment and we will suggest to place another detector of the same 
characteristics at a closer distance from the source. 

From geometrical considerations, and using the fact that the Earth is curved, a detector located on the Earth 
surface, on the vertical line of the central axis of the beam, is off-axis by a small angle, 9 mm . This is the minimum 
off-axis angle at a given distance for a given beam configuration, which for £F,N,on-axis <C R, reads 

/„ s ion — axis — -i'F.N ,„ „n 

(#min)F,N ^ , (.6.6) 

where L n-axis = 735 km is the baseline for the on-axis detector (MINOS), R is the Earth radius, and we have 
neglected terms of order (£F,N,on-axis/-R) 3 - A different (larger) off-axis angle at the same distance can be achieved 
by placing the detector slightly outside the vertical of the beam. We present here the possible locations for the near 
detector. 

It turns out that the peak energy in the CC neutrino event spectrum is well fitted by the parametrization [70J 

1900 

iW = GeV, (3.4) 

with 9 the off-axis angle in units of mrad. Then one has to solve for the detector locations that have a constant ratio 
of L/E ~ 810 km/2.3 GeV = 352 km/GeV. From Eq. I|3.4|) , it is clear that we can write 9 as a function of the baseline 
L for constant L/E, which reads 

"^"Ht)'"^)'"- 1 )- <"» 

Once the possible values of the off-axis angle 9 and the associate peak energy -E pea k are determined for a given 
distance L, it is absolutely necessary to check whether or not the geometry of the Earth and the NuMI beamline 
allow the different configurations to be a reality From Eqs. (|3.3|l and (|3.5II . the condition for this to be possible 
is 9(L) > 9 m [ n . In particular, there are no sites between 300 and 400 km which give an L/E ratio of 352 km/GeV. 
We have explored four possible detector locations: 200, 434, 500 and 600 km [7(|. Here we only present the results 
for the detector locations at 200 km and 434 km, since the distinction of the neutrino mass hierarchy is easier for 
shorter baselines, as we will show below. These configurations correspond to off-axis angles of 42.2 and 23.6 mrad, 
respectively. We have used the neutrino and antineutrino fluxes available at different off-axis configurations at the 
NO^ A far site (810 km) and obtained the fluxes at the short baselines by a simple rescaling of the fluxes at the far 
distance 6 . 

In the old proposal ^(j> the NO^A far detector at 810 km is a 50 kton tracking calorimeter, and the efficiency for its 
accepting a v e event from — ► Vf, oscillations is approximately 21%. We have explored here the possibility of using a 
50 kton liquid argon TPC detector jx^, for which the efficiency to identify v e CC interactions is 90% (i. e., basically 
perfect efficiency) and that the background is dominated by the intrinsic v e and v e components of the beam. The high 
detection efficiency of a liquid argon detector makes its statistics equivalent to that of a conventional detector, with 
the beam power upgraded with the proton driver. In addition, the physics potential of these detectors is remarkable, 
as supernova neutrino detectors, for proton decay searches 7 and for studies of neutrinoless double beta decay (see 
Ref. [721 and references therein). 

We have assumed that the number of protons on target per year is 3.7 x 10 20 Hjj (18.5 x 10 20 pot/yr with the 
Proton Driver) and five years of neutrino running. The revised NO^A proposal |4l) suggests a number of protons on 
target per year which has been upgraded to 6.5 x 10 20 ( 25 x 10 20 with a Proton Driver) and a 30 kton detector with 
24% efficiency. 



6 The short-baseline fluxes at 200 km can be easily obtained from the fluxes at 735 km and 30 km off axis l7ll . The fluxes at 434 km 
have been computed from the ones at 810 km and 18 km off-axis. 

7 If protons decay primarily into kaons (p — > K + + u), the detection efficiency in water-Cerenkov detectors is relatively low, for kaons at 
these energies are below threshold. 
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FIG. 1: (a) Approximate probability difference, Eq. {4-0$ , as a function of the neutrino energy for normal (black thick solid 
curve) and inverted hierarchies (dashed thick black curve) for sin 2 #13 = 0.02. We also depict the exact computation of the 
probability difference for the two hierarchies and different values of 5. From smaller to larger values of \T>\: 5 = ^ (blue), 
(cyan) and ^ (red) (color online). In the x-axis we specify the distances that we have explored as possible locations for the 
near detector, related to the neutrino energy by Ln = E^Lf / Ef . In this study we present the results for Ln = 200,434 km. 
(b) Same as (a) but for sin 2 2813 = 0.07. 

A. Oscillated statistics 

For a given value of the oscillation parameters, we have computed the expected number of electron events, Ni 
detected at the possible locations I = N. F (near/far sites). The observable that we exploit, A>. reads 

Ni,± = / $tj,(E v ,L) a v {E v ) P U (E„, L, fl 13 , 6, Am§ ls a) dE v (3.6) 

JB min 

where the sign +(— ) applies for the normal (inverted) hierarchies and a is the set of remaining oscillation parameters: 
$23, #12, Am|i and the matter parameter A (which depend on the baseline under consideration), which are taken to 
be known; <1>^„ denotes the neutrino flux and <7„ the cross sections. The neutrino fluxes are thus integrated over a 
narrow energy window, where -E m i n and E max refer to the lower and upper energy limits respectively. 

For our analysis, unless otherwise stated, we will use a representative jAm 2 ^ = 2.4 x 10 -3 eV 2 , which lies within 
the best-fit values for the Super-Kamiokande and K2K experiments. However, we will also show how the results 
change for different values of this parameter. For the rest of the parameters, $23, Am^ and #12, we will use the best 
fit values quoted in the introduction. We show in Appendix^ the expected number of signal and background events 
at the far location (L = 810 km) and at two of the possible near locations (L — 200 km and L = 434 km) for both 
hierarchies and four central values of the CP phase S = 0, ^, 7r, ^ and sin 2 2# 13 = 0.058. 



IV. MATTER EFFECTS AND THE TYPE OF HIERARCHY 



The first task of a long-baseline neutrino experiment should be to measure the small 6*13 angle of the neutrino 
mixing matrix. Adding one detector to the experiment would increase the statistics and then the sensitivity to this 
mixing angle. Here, we will assume that the value of #13 is within the sensitivity of next-generation long-baseline and 
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FIG. 2: fa) 27ie dashed curves depict the probability difference up to second order in Eq. \4-'^ , as a function of the neutrino 
energy for normal (lower curves) and inverted hierarchies (upper curves) for sin 2 2#i3 = 0.02 together with the exact computation 
(solid curves). From smaller to larger values of\T>\, we plot three different values of S — -| (blue), (cyan) and ^ (red) (color 
online). In the x-axis we specify the distances that we have explored as possible locations for the near detector, related to the 
neutrino energy by Ln = E^Lf / Ef . In this study we present the results for Ln = 200 and 434 km. (b) Same as (a) but for 
sin 2 2#i3 = 0.07. 



reactor neutrino experiments. We will focus on the possibility to discriminate the type of neutrino mass ordering and 
will show that, by considering the experimental setup described above, the improvement with respect to the current 
NO^A proposal is remarkable. The study of the enhanced sensitivity to the value of 9m and to the CP-violating 
phase 5 for the experimental configuration presented here will be performed elsewhere |7J|. 

In order to study matter effects, we will consider the probability of v^-v e oscillations in matter at two different 
lengths of the baselines, (near detector) and L F (far detector). And, as mentioned above, since the sensitivity to 
the mass hierarchy is optimized for L/E constant, we will consider both baselines so that we keep the same ratio L/E 
at both detectors. We compute the quantity 

V= P(L N )-P(L F ) 

~ P(L N )+P(L F ) ' ^ 

i.e. the normalized difference of the oscillation probabilities computed at the near and far locations. By using the 
approximate formula Eq. (|2.1|) and expanding up to first order in matter effects, AL <C 1 and A <C A 13 , neglecting 
the solar neutrino mixing parameters and keeping terms of 0(A9is), T> reads 

A N L N - A F L F ( 1 1 \ , s 

1 ~ 2 l(A 13 L/2) tan(A 13 L/2)J ' 1 ' > 

Let us note that the leading term in Eq. (|2.1|l is proportional to sin 2 (Ai3L/2) and cancels out in T>\ because L/E is 
the same for the near and far detector sites. It is clear from Eq. (|4.2(l that the probability difference T>\ changes sign 
if the neutrino mass spectrum is normal or inverted, since it depends on the sign of Am^. 

The effects due to the CP-phase 5, as well as those due to the solar neutrino mixing parameters, are subleading 
for large enough values of sin 2 6*13 > 0.01, region which is within the range expected to be explored by the NO^A 
experiment. In Fig. \I\we depict the neutrino energy dependence of the quantity T>i, Eq. H4.2J) . for the two possible 
hierarchies and we compare the results with the ones obtained using the exact oscillation probabilities for three 
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FIG. 3: (a) Bi-neutrino event ellipses at the short (200 km) and at the far distances (810 km) for normal (lower blue) and 
inverted (upper red) hierarchies. The experimental setup considered here is 3.7 x 10 20 protons on target per year, a 50 kton 
detector with perfect efficiencies at each location and five years of data taking. From bottom up, the ellipses correspond to 
sin 2 20i3 = 0.0003, 0.0004, 0.006, 0.008, 0.001, 0.005, 0.01, 0.02, 0.03, 0.04, 0.06, 0.07, 0.095 and 0.11. (b) Same as (a) but with the 
near detector located at 434 km. 



different values of the CP-phase S = 0, ^, ^L. We show our results for two different values of sin 2 28i 3 = 0.02, close 
to the sensitivity limit for the NO^A experiment and sin 2 29 13 — 0.07, close to the present upper bound, Eq. I|1.3|l . 

For large enough values of sin 2 26i 3 , the contribution of CP effects to T>\ is never larger than 20-30% (see Fig.^l. It 
can be shown that, at leading order, the second term of the r.h.s. of Eq. 1)2. l[l . which controls the CP-violating effects, 
depends only on L/E. Therefore, in the normalized difference T>, the leading CP-violating contribution cancels out 
and CP-violating effects can be treated as a perturbation to the dominant contribution from matter effects. 

The expression for V with terms up to 0(A A12/A13), O(A0l 3 ) and 0(A 2 ) reads 8 : 

2 x \ Ai 3 tan6» 2 3 sm20i 3 sin(Ai 3 L/2) C ° S 1 13 1 ' 2 ^ 13 



\ A\Ll-A\Ll ( 1 1 \ 

2 4 UAiaW sin 2 (A 13 L/2)J ' { ' ' 

where the correction due to the CP-phase also changes sign with the hierarchy type. For instance, for normal hierarchy, 
the correction to T>\ due to the first term in the r.h.s. of Eq. I|4.3|l will increase |X>i| if 7r/2 < 5 + \ Ai 3 \L/2 < 3ir/2 
whereas for inverted hierarchy it will decrease [D\\ if ir/2 < 5— \Ai 3 \L/2 < 3ir/2. However, the second-order correction 
due to matter effects does not depend on si<;n(Am§ 1 ) and it is always negative. 

The probability difference up to second order, Eq. 1)4. 3[l . is depicted in Fig. for the three different values of the 
CP-phase 6 — 0, ^, ^ and for the two possible hierarchies. We repeat the same exercise as in Fig. ^ and show our 
results for two different values of #13: sin 2 2#i 3 = 0.02 and sin 2 20i 3 = 0.07. 

We have thus shown that, under this special experimental configuration, determining the type of hierarchy requires 
only establishing whether T> is positive or negative, being this measurement free of other degeneracies; the corrections 
due to the rest of the neutrino mixing parameters cannot flip this sign. This implies that the determination of the 



In order to obtain the term 0(A9^ 3 ) one has to use the expression for the probability up to O(0f 3 ) (see Ref. |2fi| V and not just Eq. 12.11 . 
For large values of 813 this term is of the order of the 0(A 2 ) contribution. 
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FIG. 4: Bi-event neutrino-antineutrino ellipses at the far distances (810 km) for normal (lower blue) and inverted (upper red) 
hierarchies. The experimental setup is the one described in Fig. but considering 5 years of neutrino running and 5 years of 
antineutrino running. From bottom up, the ellipses correspond to sin 2 29n = 0.001,0.005,0.01,0.02,0.03,0.04,0.06,0.07,0.095 
and 0.11. 



type of hierarchy by using only the neutrino channel with two detectors at different baselines suffers no degeneracies 
with other parameters. In addition, the requirement of just only the neutrino channel (with two detectors, though) 
will allow the number of years of data taking (and systematic uncertainties) to be reduced. 

In order to illustrate the potential of the combination of two long-baseline detectors operating at the same L/E 
in extracting sz<7n(Am| 1 ), we show here the neutrino bi-event plots at the far and at the short distance detectors 
and we compare these results with the neutrino-antineutrino bi-event plots at just one fixed distance. In Fig. [3] we 
show the neutrino bi-event curves for the two hierarchies at the short and at the far distance, considering two short 
baselines (200 km and 434 km) and a fixed long-baseline (810 km) for different values of sin 2 26*13. As it is clearly seen 
from these plots, for the case of constant L/E at both detectors, the ellipses collapse to a line and those obtained if 
the solution is that of sign(Am'^ 1 ) > no longer overlap with those for si<?n(Am| 1 ) < 0. Notice also that the slope 
for the normal hierarchy "pencil" is smaller than that for the inverted hierarchy, because of the larger matter effect 
for larger baselines (for neutrinos). 

If instead we consider the commonly assumed configuration of using a single detector, and running first in the 
neutrino mode and then in the antineutrino mode, the ellipses overlap for a large fraction of values of the CP-phase 
S for every allowed value of sin 2 29i^. This makes the determination of sign(Am 2 1 ) extremely difficult, i. e., the 
si(7n(Am 2 1 )-extraction is not free of degeneracies. The former case is depicted in Fig. 01 where we have considered 5 
years of neutrino and antineutrino running. Therefore, even with half of the time of data taking, placing two detectors 
could resolve the type of neutrino mass hierarchy much more easily than the standard approach. 

In what follows, we will provide a detailed study of the sensitivity to the sig^Am^) in the sin 2 2613S plane. 

In order to compute the sensitivity to the mass hierarchy, as a first approach, we have constructed a measurable 
integrated asymmetry: 

a + =w? , (4 . 4) 

+ {N/N D } N + {N/N } F ' v ' 

where N is the number of v e induced events in the presence of oscillations in the normal hierarchy scheme and N a 
is the expected number of charged-current interactions in the absence of oscillations at the near (N) and far (F) 
detectors. One can compute the equivalent integrated asymmetry assuming an inverted hierarchical scenario, A-. 
If Nature has chosen, for instance, a positive value for the atmospheric splitting but the data analysis is performed 
assuming the opposite sign, the sensitivity to the sign resolution reads 
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FIG. 5: (a) Sensitivity to the sign of the atmospheric mass square difference (for lAmf^ = 2.4 x 10~ 3 eV 2 ) as defined in 
Eq. including the systematic errors induced by the uncertainties on the atmospheric mixing parameters and on the matter- 

parameter A, exploiting the data from a short-baseline off-axis detector located at 200 km and from the long-baseline off-axis 
detector at 810 km. The experimental setup assumed to obtain the solid curves is 3.7 x 10 20 protons on target per year, a 50 
kton detector with perfect efficiencies at each location and five years of data taking, whereas to get the dashed curves we have 
added to the former statistics a factor of 5 (proton driver case), (b) Same as (a) but combining the data from the long-baseline 
with a short-baseline off-axis experiment located at 434 km. 



i. e., the difference between the two asymmetries divided by its error. We have studied the case of statistical errors 
(adding backgrounds) as well the impact of the uncertainties on the remaining oscillation parameters by computing 
this systematic error on the asymmetry using the standard error propagation method. The errors on the solar 
parameters can be safely neglected. The errors considered here for the atmospheric mixing parameters Anijj and 
sin 2 2^23 are at the level of 5% and 2% respectively El], UM- For the matter parameter we take the conservative 
assumption AA/A = 5% |75| . 

We present the 95% C.L. sensitivities to the neutrino mass hierarchy resolution in Fig. [SI where Nature's solution for 
the neutrino mass spectrum has been chosen to be the normal hierarchy. We have studied both possibilities, normal 
and inverted hierarchy, as Nature's or "true" solution. We find that the conclusions do not change in a significant 
way when considering the inverted neutrino mass spectrum. 

Two possible combinations of the data at the far, "fixed" NO^A off-axis experiment have been explored: with 
the data from a near detector located at 200 km and with those from a near detector located at 434 km. All these 
experiments would exploit the NuMl neutrino beam in an off-axis mode. From the results depicted in Fig. EJ one can 
observe that the best option for the location of the second, near detector, would be 200 km: the hierarchy could be 
determined regardless of the value of the CP-phase S down to values of sin 2 2#i3 = 0.05 or sin 2 2#i3 = 0.02 with the 
Proton Driver option. For a longer baseline, 434 km, the type of hierarchy can be uniquely determined independently 
of 8 for sin 2 26>i3 = 0.06 and sin 2 20^ — 0.02, without and with the Proton Driver, resp ectively. 

For comparison, we consider the study performed in the revised NOz/A proposal |4l| on the sensitivity to the mass 
hierarchy. The analysis considers 3 years of running in the neutrino and antineutrino modes, with a 30 kton detector 
located at a baseline of 820 km and 12 km off-axis. A value of Am^ = 2.5 x 10~ 3 eV 2 was chosen. NOi/A alone 
can resolve the sign of Am^ at the 95% C.L. for only 10% of the values within the range of the CP-phase S if 
sin 2 20i3 = 0.04. Even if sin 2 26>i3 = 0.1, a value that is very close to the present upper bound, NO^A could uniquely 
determine the type of hierarchy for only 40% of the values of the CP-phase 5 ^lj ■ Even with the Proton Driver, NOi/A 
cannot resolve the sign of the atmospheric mass difference for 80% (40%) of the values of 6 if sin 2 26q 3 < 0.02(0.1) ED. 
In addition, the possibility of adding a second 50 kton detector at ~ 700 km and 30 km off-axis, is discussed El], [63 . 
In this case, the resolution of the mass hierarchy down to sin 2 2#i3 = 0.02 can be accomplished after 12 years of 
NO^ A data plus 6 years with the second detector equally split between neutrinos and antineutrinos and with Proton 
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FIG. 6: (a) Results of the \ 2 analysis to the sign of the atmospheric mass difference extraction versus sin 2 29i3, by exploiting the 
data from afar long-baseline experiment at 810 km and from a short-baseline experiment at 200 km, for |Am|i| = 2.4x 10 -3 eV 2 . 
The corresponding 90%, 95% and 99% C.L.s are shown. As a function of sin 2 2#i3, we depict the maximum (solid line) and 
minimum (dashed line) of A\ 2 , which are obtained for different values of 5 depending on sin 2 2#i3. (b) Same as (a) but with a 
Proton Driver. 



Driver. This is comparable to what the Super-NCvA setup proposed here can achieve after 5 years of running with 
only neutrinos and with Proton Driver (see Fig. EJ). However, a direct comparison needs to take into account the 
different choices of detectors, the numbers of protons on target, the value of Arn^j used in the analysis, and the fact 
that an optimization of the experiment has not been performed in our case. 

We have performed an independent \ 2 analysis of the data on the sin 2 26\^-5 plane. At a fixed value for the former 
two parameters, the \ 2 m the combination of two baselines (near and far sites) reads 

Xw = J2 - ^,±)C«'(^',± - N t >, ± ) , (4.6) 

w 

where the + (— ) sign refers to normal (inverted) hierarchy and C is the covariance matrix, which for the particular 
analysis considered in the present study contains only statistical errors. The experimental "data", A/^±, are given by 



Af £t± = (N e<± + N u ) - N be ,± , (4.7) 

where we have considered that the efficiencies are fiat in the visible energy window, Nw are the background events 
and () means a Gaussian/Poisson smearing (according to the statistics). We have assumed that nature has chosen 
a given sign for Aro 2 ^, but the data analysis is performed with the opposite sign. We show the results of our x 2 
analysis in Fig. [HI where we plot the sensitivity to sign{JS.m\ 1 ) after the combination of the data from the long-baseline 
(810 km) with the data from a short baseline at 200 km. For every value of sin 2 29\z we find the two values 9 of 5 that 
maximize and minimize A\ 2 - We depict in Fig. [5] the value of the minimum and maximum A% 2 versus sin 2 2613. If 
sin 2 20i3 > 0.04 a misidentification in the sign of the atmospheric mass difference can be excluded at 95% C.L. in the 
most pessimistic situation. There exists, however, a large number of values of 6 at which this misidentification can 
be excluded at a confidence level larger that 99%. With a proton driver, if sin 2 29^ > 0.02, our analysis shows that 
it is possible to determine sign(/S.m\ 1 ) at the level of the 99% C.L. for the full range of 8. The results from the \ 2 
analysis presented here agree with the previous study of the asymmetries. 



The values of 6 at the maximum and minimum of A\ 2 are in general different for different values of sin 2 28is . 
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FIG. 7: Sensitivity to the sign of the atmospheric mass square difference as defined in Eq. \4-5\j , for different values D/IAm^l 
= 2.0 x 10~ 3 eV 2 (dotted red line), 2.4 x 10~ 3 eV 2 (solid blue line), and 3.0 x 10~ 3 eV 2 (dashed red line). We have included 
the systematic errors induced by the uncertainties on the atmospheric mixing parameters and on the matter parameter A, and 
exploited the data from a short-baseline off-axis detector located at 200 km and from the long-baseline off-axis detector at 810 km. 
The experimental setup assumed to perform the solid curves is 3.7 x 10 20 protons on target per year, a 50 kton detector with 
perfect efficiencies at each location, and five years of data taking. 



A. Dependence on 



|Am!i| 



In the previous section, we have assumed the knowledge of lAmJ-J with a 5% uncertainty. In particular, we have 
taken lAmJjJ = (2.40 ± 0.12) x 10~ 3 eV 2 . Although the former level of precision is expected to be achieved by the 
time this experiment could turn on |40l l74| , currently the value of the atmospheric mass difference is not known with 
that level of accuracy. Thus, it is important to investigate how the results presented above change if IAto^I happens 
to be different from the previously assumed value. 

In Fig. [7| we have depicted the sensitivity to the sign of the atmospheric mass square difference determination as 
defined in Eq. (|4.5() . for three different possibilities for its absolute value, IAto^I = 2.0 x 10 _3 eV 2 (dotted red line); 
2.4 x 10 _3 eV 2 (solid blue line) and 3.0 x 10 _3 eV 2 (dashed red line). We have assumed a near off-axis detector located 
at 200 km. As can be seen from the figure, the larger the value of lAmfjJ, the better the sensitivity to the type 
of neutrino mass hierarchy. The reason can be easily understood from Eq. I|4.2|l . The asymmetry depends on the 
factor 1/x — 1/tana;, where x cx Am 2 ^ (recall that L/E is the same for both detectors). This function increases 
monotonically as x increases, and therefore the asymmetry is larger for larger lAm 2 ^, which correspondingly means 
better sensitivity. 

On the other hand, it must be remarked that since the CHOOZ 0] bound is weaker for small values of lAmJJ, 
the loss in range for is not as large as one would naively think from Fig. [7| 

In case the actual value of lAm 2 ,^ happens to be in the low side of the currently allowed range 0, a possible 
solution could be to adopt a larger L/E, which can be accomplished either by considering longer baselines or larger 
off-axis angles, i.e., smaller energies. However, both possibilities imply the reduction of the flux at the detectors, so 
a compromise must be achieved. Nevertheless, if #13 is very small and lArra 2 ^ is also small, then a longer run in the 
neutrino mode would unavoidably be needed in order to increase the statistics. All in all, a detailed analysis would 
be required to find the best possible configuration as a function of lAm^ |z3- 
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V. CONCLUSIONS 



Establishing the type of neutrino mass hierarchy — be it normal or inverted — plays a crucial role in our understand- 
ing of neutrino physics. Future long-baseline experiments will address this fundamental question. Typically the deter- 
mination of the hierarchy in the proposed experiments suffers from degeneracies with other CP-conserving and CP- 
violating parameters, namely #13, S and #23- Resolving such degeneracies in one experiment is v ery challe nging, if not 
impossible. Different strategies have been studied, e. g., by combining more than one experiment |59ll60ll6lL 163.163 . [65J. 
using more than one detector [HH l55l I5H l57l l58j . or using additional information from atmospheric neutrino data |64j . 

In the present article, we have presented a method for establishing sig^Am^), free of degeneracies, by using only 
one experiment running in the neutrino mode alone. We have considered an experimental setup with two neutrino 
detectors placed in a special off-axis configuration. It is known that off-axis spectra are well peaked at a certain 
neutrino energy, which depends on the angle from the central axis of the beam. This allows both detectors to be 
located in such a way that they have the same L/E. We have shown that very interesting synergy effects show up 
with this special configuration for which vacuum oscillation phases are the same at both sites, stressing the different 
matter effects. These are manifest when comparing the bi-event neutrino-neutrino (Fig. |3J) and bi-event neutrino- 
antincutrino (Fig. plots above, for which a clear distinction of the type of hierarchy is possible for the former 
regardless of the value of 5, but more challenging for the latter. 

We have considered a normalized difference, T> (see Eq. (|4.1|l ). between the neutrino oscillation probabilities at 
two baselines. At leading order, the sign of T> depends only on matter effects, i. e., on the type of neutrino mass 
hierarchy, while other parameters are subdominant. Although CP-violating terms can have a sizable contribution 
for small values of sin 2 2#i3, they cannot change the sign of V. This implies that the determination of the type of 
hierarchy, exploiting the method discussed here, suffers no degeneracies from other parameters. We have confirmed 
such a result by performing an analysis of the sensitivity to the mass hierarchy in a specific experimental setup, which 
we have named Super-NO^A. We propose to use the NuMI beam in the neutrino mode and two detectors, one at the 
far distance proposed by the NO^A Collaboration, Lp ~ 800 km, and the other with a shorter baseline, Ln ~ 200 km 
(434 km), with the energy tunned to E N = E F L N /L F . The selection of the short baseline must be done in such a way 
that it is possible to place a detector at the precise off-axis angle, in order to get that particular energy at the peak 
of the spectrum. Because of the Earth curvature, the near detector, located on the Earth surface and on the vertical 
of the on-axis beam, is off-axis by a small angle, which is the minimum possible off-axis angle at that distance. This 
implies that not all different configurations, such that L/E is the same at both sites, are possible. In particular, there 
are no sites between 300 and 400 km which give an L/E ratio of 352 km/GeV. 

By considering the integrated asymmetry, Eq. (|4.4|) , we have shown that a suitable baseline for the second detector 
to determine the type of neutrino mass hierarchy corresponds to L^ ~ 200 km, which enhances matter effects without 
the need of too low energies to keep the same L/E at both detectors. We have shown in Figs. [S] and El that this can be 
achieved at 95% C.L., regardless of the value of S, for sin 2 2#i3 > 0.05 for a conventional beam and for sin 2 29is > 0.02 
with a Proton Driver. Similar results can be obtained for a slightly longer baseline, e. g., 434 km. We have also 
performed an independent \ 2 analysis of the simulated data on the sin 2 2613-8 plane, which confirms our previous 
results. This is in contrast with the sensitivity of the proposed NOz^A experiment. At the 95% C.L., only for 10% of 
the values of the CP-phase 5, NOz/A can resolve the type of neutrino mass hierarchy if sin 2 2#i3 = 0.04, considering 
three years of neutrino plus three years of antineutrino running |40l l4l| . 

In Fig. we have also shown the sensitivity to the type of neutrino mass hierarchy for three values of lAm 2 -^ for 
the adopted configuration, and as can be seen from it, if lAm^l lies in the low side of the presently allowed range 0, 
another configuration [73j or more statistics might be needed. 

For our simulations we have considered two 50 kton liquid argon detectors. In addition to the off-axis experiment 
detailed here (and in general, to any long-baseline neutrino experiment), the physics potential of liquid argon detectors 
is remarkable. They can also be used as supernova neutrino detectors, for proton-decay searches, and for studies of 
neutrinoless double beta decay (see Ref. |72| and references therein). 

It is important to notice that the use of the neutrino mode alone would allow a reduction of the number of years 
of data taking if compared with the standard approach of running in the neutrino and then antineutrino modes. In 
addition, having two identical detectors and only one beam reduces the systematic uncertainties. 

Thus, we have presented an improved off-axis experiment with respect to the proposed NOz^A experiment with a high 
sensitivity to the type of neutrino mass hierarchy (free of degeneracies) even with only a neutrino run. The improved 
capabilities for measuring the value of #13 and the CP-violating phase, as well as other possible configurations, will 
be studied elsewhere 1731. 
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APPENDIX A 



We present in this appendix the computed charged-current neutrino event rates for the NuMI beam and different 
locations of a 50 kton liquid argon detector. As a matter of illustration, we show these event rates in Tables U ITT1 
and IIIII for a given value of sin 2 2#i3 = 0.058, a given value of lAmU = 2.4 x 10~ 3 eV 2 , and for four different values 
of the CP-phase, 5 = 0, 7r/2, it, 3tt/2. In these tables we show the unoscillated ^-like events, the expected oscillated 
z; e -like signal and the v e intrinsic background, which mainly comes from /i decays. 



sin 2 20i3 


6 


Ami, (eV) 


Vf± (unoscillated) 


v e (signal) 


v e (intrinsic background) 


0.058 





0.0024 


16322 


481 


59 


0.058 


7r/2 


0.0024 


16322 


526 


59 


0.058 


7T 


0.0024 


16322 


367 


59 


0.058 


3tt/2 


0.0024 


16322 


323 


59 


0.058 





-0.0024 


16322 


264 


59 


0.058 


7r/2 


-0.0024 


16322 


398 


59 


0.058 


7T 


-0.0024 


16322 


353 


59 


0.058 


3tt/2 


-0.0024 


16322 


219 


59 



TABLE I: Calculated charged-current neutrino event rates (signal and backgrounds) for NOvA (baseline of 810 km, 10 km off- 
axis), for 3.7 x 10 20 pot/yr, 5 years running at a 50 kton far detector. We have computed them for |Ato§]J = 2.4 x 10 -3 eV 2 , 
sin 2 20i3 = 0.058 and four values of 5 = 0, ?,7r, 4^. The remaining oscillation parameters are AmJi = 8.0 x 10~ 5 eV 2 , 
sin 2 0i2 = 0.31, sin 2 2023 = 1 and the matter parameter A = y/2Gpn e = 1.064 x 10 -4 eV 2 /GeV. The energy window is [1.8,2.8] 
GeV. The neutrino spectrum peaks at 2.3 GeV. 



sin 2 20i3 


S 


Ami, (cV 2 ) 


z/ M (unoscillated) 


u e (signal) 


v e (intrinsic background) 


0.058 





0.0024 


13285 


342 


81 


0.058 


tt/2 


0.0024 


13285 


373 


81 


0.058 


7T 


0.0024 


13285 


257 


81 


0.058 


37r/2 


0.0024 


13285 


226 


81 


0.058 





-0.0024 


13285 


228 


81 


0.058 


7r/2 


-0.0024 


13285 


335 


81 


0.058 


7T 


-0.0024 


13285 


301 


81 


0.058 


3tt/2 


-0.0024 


13285 


193 


81 



TABLE II: Same as Table\]\ for a baseline of 434 km (10.23 km off-axis). The matter parameter A = y/2Gpn e = 8.93 x 
10~ 5 eV 2 /GeV, and the energy window is [0.8,1.8] GeV. The neutrino spectrum peaks at 1.3 GeV. 
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sin 2 20i3 


S 


Ami, (eV) 


i/ M (unoscillated) 


v e (signal) 


v e (intrinsic background) 


0.058 





0.0024 


8253 


163 


83 


0.058 


tt/2 


0.0024 


8253 


177 


83 


0.058 


7T 


0.0024 


8253 


122 


83 


0.058 


37r/2 


0.0024 


8253 


108 


83 


0.058 





-0.0024 


8253 


123 


83 


0.058 


7r/2 


-0.0024 


8253 


177 


83 


0.058 


7T 


-0.0024 


8253 


160 


83 


0.058 


3tt/2 


-0.0024 


8253 


106 


83 



TABLE III: Same as Tabled for a baseline of 200 km (8.44 km off-axis). The matter parameter A = y/2Gpn e = 3.83 x 
10 -5 eV 2 /GeV, and the energy window is [0.2,1.2] GeV. The neutrino spectrum peaks at 0.7 GeV. 
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